Decoherence and the quantum-classical limit in the presence of chaos. 
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We investigate how decoherence affects the short-time separation between quantum and ciassicai 
dynamics for ciassicaily chaotic systems, within the framework of a specific model. For a wide 
range of parameters, the distance between the corresponding phase-space distributions depends on 
a single parameter \ that relates an effective Planck constant h e g, the Lyapunov coeffficient, and the 
diffusion constant. This distance peaks at a time that depends logarithmically on h B ff, in agreement 
with previous estimations of the separation time for Hamiltonian systems. However, for \ i$ F the 
separation remains small, going down with ftjiff > so the concept of separation time loses its meaning. 
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One of the most subtle problems of quantum mechanics 
is the description of the classical world, particularly for 
classically chaotic systems. Even for initial states that 
are classically allowed, one expects that the dynamics 
of the quantum and the corresponding classical system 
should differ, after some time. Indeed, while for linear 
systems the Wigner distribution can be shown to obey 
the same dynamical equation than the classical phase- 
space distribution, non-linearities will eventually set the 
two distributions apart. 

For a classically chaotic system, the separation time 
can be very short, due to the exponential stretching 
of the distribution, for positive Fyapunov coefficients, 
which quickly allows the distribution to explore the non- 
linearities of the system, even if the linear dimensions of 
the initial wavepacket are smaller than the typical non- 
linear scale of the problem. The separation time in this 
case has been shown by many authors 0, @ to scale as 
ln(l/?ioff), where h s — Ti/S, and S is a typical action 
of the system. Thus, even when h e s <SC 1, the separa- 
tion time can still be small, as compared to a typical 
evolution time of the system. Such logarithmic law may 
pose a problem to the quantum-classical correspondence 
of macroscopic objects, leading to consequences that con- 
tradict observation [3|. Reconciliation of theory and ob- 
servation is provided by the irreversible coupling of the 
system with a reservoir, which leads to the elimination of 
quantum signatures. This has been investigated numeri- 
cally and analytically for several models U 0, @ • 

In the presence of the environment, one expects that 
the logarithmic law should not hold anymore. In Ref. 
it was shown that, for the kicked harmonic oscillator, dif- 
fusion helps to decrease the difference between quantum 
and classical variances, so that, if one defines any arbi- 
trary value of this difference as the "critical separation 
value" , the separation time becomes infinite for a suffi- 
ciently large diffusion coefficient, that is, the difference 
remains always smaller than the critical value. In that 
work, however, the dependence on the relevant parame- 
ters of the difference between classical and quantum dy- 



namics could not be elucidated. 

In this paper, we derive for the kicked harmonic os- 
cillator 0, the precise dependence of the separation 
between quantum and classical distributions on the pa- 
rameters that control macroscopicity, noise, and chaotic 
behavior. These factors, for a wide range of their val- 
ues, can be combined in a single parameter, which will 
be shown to govern the quantum-classical transition, as 
conjectured in Ref. [5j for general chaotic systems. 

The relevant dimensionless parameters correspond to 
the diffusion coefficient D, the kicking strength K, and 
an effective Planck's constant ft e ff- We show that, in the 
chaotic regime, for any finite diffusion coefficient, and in 
the semiclassical limit h c g <C 1 , the distance between the 
two distributions, defined as the integral of the magni- 
tude of their difference over all phase space, is propor- 
tional to x — Kfi-cfi /4-D 3 / 2 , as long as x < 1. In this 
regime, the time for which this distance peaks is shown 
to be a logarithmic function of h e g, as estimated previ- 
ously for Hamiltonian chaotic systems. However, in this 
case the concept of separation time is not meaningful any- 
more, since the two distributions remain close together 
throughout the evolution. 

The Hamiltonian of the kicked harmonic oscillator is 
defined as: 

pi i 00 
H = — + -mv 2 Q 2 + A cos{kQ)^2s(t-nT). (1) 

This Hamiltonian can be shown to describe the center- 
of-mass dynamics of an ion in a one-dimensional trap 
submitted to a sequence of standing-wave laser-pulses, 
off-resonance with a transition between the electronic 
ground state and another internal state ||- The wave 
number k in Eq. is the projection along the trap axis 
of the corresponding (identical) wave vectors of two op- 
posite propagating pulses with oblique incidence. The 
high degree of control in ion experiments [10|], plus the 
possibility of engineering several kinds of reservoir for the 
center-of-mass motion greatly stimulates the inter- 
est in using this system for testing fudamental features 
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of the quantum-classical transition. The corresponding 
phase space is unbounded, and no long-time localization 
occurs, as opposed to the kicked rotator, which has been 
subject to experimental test concerning the dependence 
of localization on noise . 

It is convenient to work with the dimensionless quan- 
tities q — kQ, p — kP/mv, and K — k 2 A/mv, so that 
[q,p] = 2irj 2 = ih c g, with r\ — kAQo — k^Jfi/2mv, and 
AQo being the width of the ground state of the har- 
monic oscillator. The dimensionless parameter r\ is the 
so-called Lamb-Dicke parameter 10f , which measures the 
ratio between the ground state width and the wavelenght 
A = 2n/k that sets the scale of the non-linearity of the 
Hamiltonian. It is important to note that, in experiments 
with trapped ions, the classical limit r\ — > can be ap- 
proximated simply by changing the angle of incidence of 
the incoming pulses, or by increasing the trap frequency. 

The time-dependent classical evolution can be de- 
scribed as a composition of a discrete map correspond- 
ing to the kick, plus a rotation in phase space: x n+ i = 
Ro K(x„) where x„ = (q n ,p n ) are the coordinates be- 
fore kick n (note that our first kick corresponds to n = 0). 
The operation K is defined as 

Qn = In , Pn = Pn + K Sm(q n ) , (2) 

and the phase space rotation R is given by: 



may be described by the Fokker-Planck equation for the 
Wigner function 
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where T is the dissipation rate. For the complete evolu- 
tion we must add to the right-hand side (r.h.s.) of Eq. (J3J) 
the unitary evolution given by Eq. In the classical 
case, the complete evolution is obtained by adding the 
corresponding Liouville term to the r.h.s. of Eq. (J5J. 

In the low-temperature limit tt, <gc 1 , the diffusion term 
becomes negligible in the semiclassical regime 77 <gc 1 , and 
one gets then purely dissipative dynamics. On the other 
hand, when n — > 00, T — > 0, with nT constant, one gets 
a purely diffusive dynamics. In this limit, Eq. JSJ) is also 
the diffusion equation for a classical distribution W cZ (x), 
as long as nTr] 2 is identified with a classical diffusion 
constant T. 

The first limit has already been investigated in Refs. 0, 
ll2T |. where it was shown that pure dissipation, in the 
sense described above, does not change the logarithmic 
law for the separation time. In fact, they show that the 
separation time ts, for an initial wave packet centered 



q n+ i = cos(jat) g+ + sin(^r)p+ , 
Pn+i = - sin(^r) q+ + cos(^r) p+ 



(3) 



For definiteness, we consider here r = T/6 (T = 27r/^ 
is the period of the harmonic oscillator) wich leads to 
the so-called "stochastic web" in phase space Q, i.e. 
a pattern of groups of stability islands with hexagonal 
symmetry immerse in a sea of chaotic trajectories. The 
strong-chaos regime corresponds to K 3> 1. 

In order to study the separation between quantum and 
classical dynamics, we shall use the quantity 



V n = J dx|VF„(x) 



W n ci (x)| 



(4) 



where W n and are respectively the Wigner and the 
classical distributions immediately before kick n. The 
normalization of the distributions is taken equal to one. 
A typical example of the variation of D n as a function 
of the number of kicks, in the absence of decoherence, 
is shown Fig. ^ Separation becomes evident as T> n gets 
larger than one. The first peak of T> n corresponds to the 
first folding of the distributions, as their widths become 
of the order of one and the system starts exploring the 
nonlinearity of the Hamiltonian. 

We discuss now how decoherence affects the behavior 
of T> n . The effect of a thermal reservoir with average 
population n, in the Markovian and weak-coupling limit, 
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FIG. 1: (Color online) Separation T> n between quantum and 
classical distributions as function of the number of kicks, for 
77 = 0.3, K — 2, for unitary evolution. The initial distribution 
is a coherent state centered around the origin of phase space 
with width Aq = Ap = 77 much smaller than the nonlinear- 
ity scale (77 1). The peak correponds to the first folding 
of the distribution, leading to the appearance of interference 
fringes in the Wigner function (upper inset). The lower inset 
displays the classical distribution. The folding appears when 
the distribution starts probing the nonlinearity. 
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around the origin of phase space, is given by 



ts^jHl/r]), 



0.015 



(6) 



where A is the logarithm of the expansion eigenvalue of 
the linearized map at the origin. With dissipation, the 
expansion eigenvalue includes a factor exp(— Tt/2), so 
in this case A = A — Tt/2, where Ao is the expansion 
eigenvalue without dissipation. If instead one takes an 
average over initial conditions in the chaotic region, then 
one can show that A is replaced by the Lyapunov coeffi- 
cient A. With dissipation, one has A = Ao — Tt/2. The 
limit Tt/2 = Ao is never attained in the chaotic regime, 
since before that the chaotic behavior disappears, being 
replaced by simple attractors in phase space @. This 
implies that dissipation by itself is not an efficient mech- 
anism to increase the quantum-classical separation time, 
as defined by Eq. |0. In fact, the logarithmic dependence 
remains, in spite of the dissipation. For the system con- 
sidered here, in the strong chaos and weak dissipation 
regime, Ao = In [(K/2) sin(^r)] and Ao = In [K sm(vT)}. 

We show now that diffusion can drastically change this 
scenario. Without a reservoir the unitary evolution of the 
Wigner function can be written as 



W„+i(x) = J dx' L(xV) W n (x') , 



(7) 



where L(x , x') is the propagator 



2Tvrf 



S{q R -q') , (8) 



corresponding to one kick plus a harmonic evolution, with 
x R = [q R (x) , p R (x)] = R _1 (x) being the phase space 
coordinates rotated with the inverse transformation of 
Eq. J2J. The Liouville evolution of the classical distribu- 
tion W^'(x) can also be written in the form of Eq. 
with the classical propagator 

L cl (x R , x') = 5 [jf - P R + K sin(g')] 6 (q R - q') • (9) 

In the classical limit 77 — s- 0, stationary-phase techniques 
[l3j | guarantee that this classical propagator is formally 
recovered from the quantum one. 

Diffusion leads to a smoothing of the propagator, 
for both classical and quantum cases. The classical 
smoothed propagator becomes 



AttD 



(10) 



where y — (p' — p R + AT sing') /2\f~D and x = 
(q> - q R ) /2VD. 

Note that when D -t Owe recover Eq. (0 . In the quan- 
tum case diffusion leads to a factor exp [—D^ 2 /rf^ in the 
integrand of Eq. JHJ . When the width of this Gaussian is 




FIG. 2: Renormalized distances between quantum and clas- 
sical distributions, with diffusion. The seven different curves 
correspond sequentially to r\ =0.1 (dashed line), 0.04, 0.02, 
0.015, 0.007, 0.005, 0.003 (dotted line) and D = 5.13 x 10" 2 
(dashed line), 4.5 x 10" 3 , 7 x 10" 4 , 3.25 x 10" 4 , 4.5 x 10" 5 , 
1.74 x 10~ 5 , 4.5 x 10~ 6 (dotted line). In all the cases, 
X = 0.017. n pca k(ry)r is the time at which the first peak of 
the distance corresponding to a given value of rj occurs. Fhe 
second peak happens at the second folding of the distribution. 



small, r) 2 /vD <C 1, the /i's that effectively contribute to 
the integration are those close to the origin. This allows 
us to use sin(/i) ps fi — /i 3 /6 in the phase of the integrand. 
Moreover, if x = Krj 4 /D 3 / 2 <C 1, the term with /j, 3 in the 
phase is small, so we can use e l ^ e+s ^ ~ e w + iS e l8 and 
then perform the /i-integration. One should note that, 
when -q < \f~K^ which is always true when K > 1 and 
i] < 1, then x < 1 implies that D > ?/ 4 . We get then 
the following approximation to the smoothed quantum 
propagator: 

Z(xVO«£ ci (xVO[l + Xsin(qO/(y)] , (11) 

where / (y) = 1/4 (y-2y 3 /3). 

The correction to the classical propagator depends only 
on x- Since |/ (y)exp (— y 2 ) | < 0.081, Eq. l(l"l"|) is valid 
under the less restrictive condition % < 1. One should 
note that x depends on the Lyapunov coefficient through 
K. For K 3> 1, the Lyapunov coefficient Ao depends 
logarithmically on K , and therefore x is proportional to 
exp(Ao), which differs from the general expression con- 
jectured in Ref. 

The fact that x rules the corrections to the classical 
propagator implies that the separation between the clas- 
sical and quantum distributions, defined by Eq. I@J, is 
also scaled by this parameter. This is shown in Fig. 
which displays the scaled separation, for a wide range 
of values of r\ and D. The horizontal axis is scaled by 
^peak(??), the number of kicks for which the first peak 
of T> n (rj) is attained. All the curves fit in the same 
scale, which establishes that, for D > (Kr] 4 ) 2 ^ 3 > rj 4 , 
which is always attainable in the semiclassical limit, in- 
dependently of the value of D, the separation between 
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FIG. 3: 

kick n - 
4.5 x 10" 
D = 4.5 



(Color online) Distributions immediately before the 
= 20, for x = 0.017. Classical: (a) r? = 0.04, D = 
" 3 , (b) n = 0.007, D = 4.5x 10~ 5 and (c) n = 0.003, 
x 10~ 6 . Wigner: (d) same parameters as in (c). 



the quantum and the classical distributions goes down 
with 77 4 , and therefore may become arbitrarily small. 
Fig. [21 also shows that the position of the peaks as a 
function of In 1/ 77 is fitted by the straight line n poa k = 
1.45 In (1/77) + 0.54, in excellent agreement with Eq. JHJ, 
since in this case 1/Aq = 1.47. Therefore, the time when 
the first peak in T> n occurs still behaves logarithmically 
with 77. However, in view of the smallness of T> n , the 
concept of separation time is not meaningful anymore. 

Note that the condition D >rf prevents us from tak- 
ing the limit D — > for every fixed value of 77. This 
restriction does not necessarily imply a diffusion strong 
enough to substantially change the chaotic classical dy- 
namics, since we are interested in the semiclassical limit 
rj 2 <C 1 . Figure El displays the classical and quantum 
phase-space portraits for \ — 0.017 and K = 2, for 
different values of D and 77. Since x IS fixed, as rj de- 
creases, so does D, implying that it is possible, even for 
D > (Krj 4 ) 2 / 3 , to have a non-trivial dynamics, which is 
not dominated by diffusion. 

Up to now, only the case x ~ 1 has been considered. 
The behavior of the distance for other values of x is dis- 
played in Fig. 01 which shows that the proportionality to 
X is valid up to V n w 1, when the concept of separation 
time becomes appropriate. 

In conclusion, we have shown that, in the macro- 
scopic limit, and in the chaotic regime, the distance 
between quantum and classical distributions scales as 
X = Krj 4 / D 3 / 2 , when x ^5 1 5 which implies that, no mat- 
ter how small D is, this distance can be made as small as 
one wants, by decreasing the effective Planck constant 77. 

We have also shown that the maximum distance be- 




FIG. 4: Maximum distance between quantum and classical 
distributions as a function of x for K = 2, and D — 4.5 x 10 -5 
(circles) or D — 4.5 x 10 -4 (crosses). The linear behavior 
(dashed line), characterized by the unit value of the slope of 
the linear region, is verified up to values of x well beyond one, 
for which the distributions are clearly separated. 



tween the two distributions is attained at a time that 
scales logarithmically with ln(l/?y). Without diffusion, 
this leads to the well-known logarithmic dependence of 
the separation time. With diffusion, and x T although 
this behavior of the peaks is still present, its use to define 
a separation time becomes meaningless. 
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